Re-reheating, late entropy injection and constraints from baryogenesis scenarios 
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Many theories of particle physics beyond the Standard Model predict long-lived fields that may 
have dominated the Universe at early times and then decayed. Their decay, which injects entropy in 
the thermal bath, is responsible for a second reheating, dubbed re-reheating, that could substantially 
dilute the matter-antimatter asymmetry created before. In this paper we analyze such late re- 
reheating and entropy dilution. It turns out that in some cases the usual analytic calculation badly 
fails if it is not rectified by some corrective factors that we provide. We also determine the parameter 
space where the entropy dilution compromises models of baryogenesis. This region can be obtained 
by imposing some generic constraints that are applicable to any baryogenesis mechanism and long- 
lived field satisfying a few assumptions. For instance, by applying them to MSSM electroweak 
baryogenesis, thermal non-resonant leptogenesis and thermal resonant leptogenesis, we obtain that 
the initial abundances of long-lived fields with lifetime longer than respectively 5 x 10 13 , 10 -2 and 
10 15 GeV -1 are strongly constrained. Similarly, the same baryogenesis scenarios are incompatible 
with large oscillations of moduli with mass smaller than 0(1O 8 ), O(10 13 ) and O(10 7 ) GeV that are 
naturally coupled to the visible sector via gravitational dimension- five operators. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 

The history of the Universe before Big Bang Nucle- 
osynthesis (BBN) is an open issue. The standard cosmo- 
logical scenario considers epochs prior to BBN to be radi- 
ation dominated. However, this assumption is question- 
able since experiments cannot put severe constraints on 
the Universe at temperature T ^> 0(1 MeV) and more- 
over several theoretical frameworks predict a modifica- 
tion of the standard cosmological picture. For instance, 
models having long lived fields that are not in thermal 
equilibrium may induce an era where the radiation energy 
is subdominant. Potentially, this epoch might emerge in 
presence of flat directions Q-balls Q, gravitinos Q, 
axinos moduli Q, which we will generically refer to 
as X field in the following. 

On the other hand, the history of the Universe at 
T < 0(1 MeV) is well established: any primordial X- 
dominated epoch must end before the onset of BBN in 
order not to jeopardize the predictions of the primordial 
element abundances [f| @|. The required return to the 
radiation dominated era may occur via X decay. This 
process, which we dub re-reheating to distinguish it from 
the first reheating happened in the inflationary epoch, 
dumps entropy into the thermal bath and may consider- 
ably dilute the pre-existing species. 

Various experimental observations can be explained by 
re-reheating. Among many we mention the possibility of 
tuning the V-decay rate to generate the right amount 
of non-thermal cold dark matter required for structure 
formation @. Moreover, the X decay can be invoked to 
wash out unwanted relics of early stages of the Universe 
(e.g., monopoles and domain walls), to circumvent the 
gravitino problem Q or to reduce an overabundance of 
thermal dark matter candidates [Icj . 



All these possibilities highlight that future detection 
of new particles might provide a link between (non- 
standard) cosmology and particle physics that might be 
misunderstood if non classical cosmological scenarios are 
not taken into account. Thus, in order to infer cosmo- 
logical issues from particle experiments, we will need to 
clarify as much as possible the picture of the Universe 
before BBN. In particular, bounds on the possible X- 
dominated epoch will be useful 0- With this goal in mind, 
in the present paper we link the puzzle of the observed 
Baryon (B) Asymmetry of the Universe (BAU) with the 
X decay. In this manner we obtain some bounds on mod- 
els involving baryogenesis and late-time entropy injection 
where the latter does not induce B violations. 

The experimental measure of the BAU is achieved by 
the BBN [Hi and CMB [H| analyses. Under the assump- 
tion that only photons and neutrinos are relativistic at 
T< 0(1 MeV), these analyses imply 



7.2 <Y e B xp x 10 u < 9.2 
8.4 <Yl xv x 10 n < 9.2 



BBN at 95% C.L. , (1) 
CMB at 95% C.L. , (2) 



where Y^ xp = N B /S* and S* (N g ) is the total entropy 
(total baryon minus antibaryon number) during the BBN 
epoch, i.e., T < Tbbn = 4 MeV @. These measures 
constrain the late-time evolution of Yb (t) = jVe(i)/ S(t) 
and, once Ns(t) is known, some bounds on S(t) can be 
inferred. In other words, by knowing when and in what 



1 Some constraints on X are deduced by noticing that the X decay 
might have left a trace in the Cosmic Microwave Background 
(CMB) HH]. Unluckily, such a signature cannot be disentangled 
from the inflaton one so that the derived bounds on X depend 
strongly on the assumed inflationary model. 
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abundance the B asymmetry was produced (i.e., by as- 
suming a given baryogenesis mechanism), one can con- 
strain the X-decay entropy injection that is compatible 
with the above BBN and WMAP measures. Alterna- 
tively, by assuming some characteristics of the X field, 
one can impose extra requirements on the mechanism re- 
sponsible for the BAU. 

Since we perform our analysis without considering any 
specific baryogenesis mechanism or any particular X 
field, our derived bounds are easily applicable to a wide 
class of models. For this reason, the generic results we 
obtain can be used as tools by which the reader can eas- 
ily estimate the parameter region where her/his favorite 
baryogenesis mechanism and entropy injecting field (ful- 
filling some requirements explained in the text) are com- 
patible. 

As illustrative applications, we show the implications 
of our results for some specific baryogenesis mechanisms. 
More precisely, we obtain the parameter region of a 
generic X field where the late-time entropy injection is 
compatible with electroweak baryogenesis and resonant 
or non-resonant thermal leptogenesis embedded in the 
Minimal Supersymmetric Standard Model (MSSM) . Sub- 
sequently we consider the specific case of X fields being 
moduli. 

The paper is organized as follows. In Section II we 
quantify analytically the entropy injection produced by 
the X decay and compare it with the numerical predic- 
tion. Section III is devoted to demonstrate how the differ- 
ent baryogenesis mechanisms constrain the lifetime and 
energy density of a generic X field and viceversa. In Sec- 
tion IV we apply our results to concrete cases of baryo- 
genesis mechanisms and long-lived X fields. Finally, we 
leave Section V for conclusions and Appendix for techni- 
cal details. 



A. Setup and numerical evolution 

The evolution of the exponentially decaying X fields 
with energy density scaling as R(t)~ 3 is described by the 
equation [lj] 



Its solution is 



px + 3Hp x = -px/r 



R 



(3) 



Px(R) = Px(Ri) ^ exp(-i/r) , (4) 

where Ri is the scale factor at some initial time ti <C t 
and H is the Hubble constant given by 



H 2 = (R/ R) 2 = (p x +p R )/(3m 2 P ) 



(5) 



with the reduced Planck mass nip = 2.4 x 10 18 GeV. 
Assuming all the X decay products to thermalize suffi- 
ciently fast implies that the total entropy S evolves as 



s={ iss - ) R px/T 



(6) 



where g(t) is the number of relativistic degrees of freedom 
in the thermal bath. In turn, S is linked to pr by the 
expression 



PR = - A 



45 ff 4 
2n 2 g(t) 



1/3 



R- 



(7) 



Therefore, the entropy dilution caused by the X decay, 



A = 



S(t > r) 
S(U) 



(8) 



II. QUANTIFYING ENTROPY INJECTION 



Let us assume the existence of a weakly-coupled field 
X with lifetime t and energy density px redshifting as 
R(t)~ 3 , where R(t) is the expansion scale factor of the 
Freedman-Robertson- Walker Universe. Before the X de- 
cay, which we assume to be mostly into components of the 
thermal bath, the radiation energy density pp falls faster 
then px by a factor R(t) roughly. Hence, if the X de- 
cay is suppressed enough, a A-dominated epoch emerges 
even though initially px is subdominant to pr. This de- 
viation from the standard history of the Universe ends 
once the X fields decay away, what may release a sizable 
amount of entropy into the thermal bath. In this section 
we quantify this entropy injection and the temperature 
at which this injection ends. 



can be obtained by solving numerically the coupled 
Eqs. (HJ)-© and imposing the initial conditions px(Ri) = 
p x and p R (Ri) = p R . 

In the numerical solution we implement g(T) assum- 
ing that the QCD phase transition takes place at T rj 
0.2 GeV Moreover, the time dependence in g(t) re- 
quired in Eq. ([5]) is derived iteratively from g(T) by solv- 
ing the above differential equations and determining T(t) 
from S(t) = 2n 2 g(T)T 3 R 3 /45. 

As an example, in Fig. 1 (dashed lines) the numerical 
evolution of pn(t), px(t) and S(t) in the Standard Model 
(SM) are shown for initial conditions Pr/px = 10 10 , Ti = 
10 9 GeV and r = 10 22 GeV -1 . The corresponding dilu- 
tion factor turns out to be A ~ 7.8. 



2 It is assumed that the X field has a number density much larger 
than its equilibrium value, or to be an unstable field, when 
3Ht < 1. 
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FIG. 1. Evolution of the energy densities pR and px and 
the total entropy S as function of time for initial conditions 
Pr/Px = 10 10 ,Ti = 10 9 GeV and r = 10 22 GeV -1 . The re- 
sult of the numerical (analytical) solution is shown by dashed 
(solid) lines. The X-dominated epoch starts at t e and finishes 
at t d . 



B. Analytic approximation 

A straightforward analytic approximation to estimate 
the entropy dilution A is based on simplifying the previ- 
ous exponential decay to an instantaneous one occurring 
at the time U ~ r [ll Fl This approximation implies 
S(t) = S{ for t < td and it yields an evolution of the 
Universe as the one described by solid lines in Fig. 1, as 
we prove now. 

Let us assume the Universe to be radiation dominated 
at the initial time ti 0- For large enough r there exists 
an equilibrium time t e (3> U) after which pn becomes 
smaller than px- Using Si = S e (quantities with the 



index i, e, d, rh are evaluated at t 
tively) one gets 



ti, te, t ^ j j respec 



,T„ 3 



Pr = Px 



fJ 



T 3 



Px 



T e = T lPx /p R . (9) 



By taking the total energy and the scale factor R to 
be constant at the moment of the decay, it turns out to 
be 



P d x =Px(t-;) = p% = p R (t+), 

and subsequently 



rh 

Pr 



9eTl 



Pr 



T 



9rh 



(10) 



(11) 



which implies an (apparent) instantaneous growth of T0. 

The time td, at which the sudden decay produces the 
discontinuity, is related to r by the definition 



(12) 



where the chosen prefactor 3/2 is motivated by the 
matter-like dominated epoch prior to the decay. Corre- 
spondingly, at t — td the X decay increases the entropy 
by the factor 



A = 



S,-l 



9rh.T? h 



Sd gdTl 

and, since Eqs. (jTUJ) and (fT2"j) imply 

40 ml 



(13) 



rri4 

1 rh ~ 



2 9 

K 9ThT 



one concludes 



j 1/4 

A = 1.38 x 10- 9 ^ ^ 



GcV" 



1/2 



(14) 



(15) 



In the rest of the paper we will use Eqs. (|T4|) and (jl5[) 
as analytic estimates of the re-reheating temperature and 
entropy dilution. 

Notice that Eq. (fT3"|) and thus (fT5")) are strongly sensi- 
tive to the evaluation of T e which derives from the equal- 
ity Si = S e that is valid only for t e <T. Hence for Awl, 
i.e. T e ~ T r h, Eq. (| 1 5|) is unreliable. On the other hand, 
Eq. (|15p turns out to be realistic for A ^ 1, as Fig. 1 
suggests and as we will show in the next section. 

Fig. 1 also highlights that the sudden decay approxima- 
tion is not useful to describe quantities at t < r at which 
the discrepancy with the numerical quantities is huge. 
In particular, it is well known that the discontinuity of T 
at t = td is an artifact of the crude approximation: ac- 
tually the temperature of the Universe always decreases 
smoothly during the decay [l5[ , as the numerical analysis 
in Fig. 1 highlights. 

Finally, notice that in order to compare the evolutions 
of the analytical and numerical quantities in Fig. 1, the 
different scale factors R(t) relative to each one of the two 
approaches are introduced. This allows to observe that 
in the numerical evolution px returns subdominant near 
t = t. what is reproduced in the analytic approximation 
thanks to the factor 3/2 chosen in the definition (fT2")) . If 
this factor were lowered, the instant of the temperature 
discontinuity would be anticipated. 



3 The present fully-analytic approach is the most straightforward 
for practical purposes, so we will compare our numerical results 
with it. For some semi-analytic approximations see Ref. Il 5ll . 

4 In practice other reliable situations can be led to this case by 
shifting and redefining the time. 



Alternatively, the instantaneous decay could be treated by main- 
taining fixed T and increasing instantaneously R. This would not 
lead to different conclusions. 
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C. Entropy production and re-reheating 
temperature 

Besides the entropy dilution, the estimates of the re- 
reheating temperature T r h obtained by the analytic and 
numerical approaches are important. However, while in 
the sudden decay approximation T r h is well defined - it 
is the temperature just after the decay - in the numerical 
method it is not. In order to extend the concept of re- 
reheating temperature to the latter case we introduce the 
convention considering that the entropy injection ends 
when its remaining variation £ is tiny. In this way in 
the numerical analysis we can define implicitly T r h as 
S(T -C T r h) I S (Trh) = 1 + £ and, after an opportune 
choice of £, we can use the bound 



T r h > Tj 



BBN 



(16) 



to avoid alterations of BBN predictions. 

To determine £ we observe that 1% of entropy in- 
jection after Tbbn roughly corresponds to the resolu- 
tion of the measurements (p} and ([2]). This means 
that even in the pathological parameter scenario yield- 
ing Y B (T = T B bn) ^ 8.4 x lCT 11 [the lower bound al- 
lowed by ([1]) and ©], a further 1% dilution does not 
lower appreciably Yg which then remains within the ex- 
perimental constraint at any T < Tbbn- For this reason, 
the requirement of not jeopardizing the BBN predictions 
can be safely replaced by the bound (fTB"|) with £ = 0.01 
(see Appendix for considerations about other choices of 
£) . Consequently, the definition of reheating temperature 
that we will use in the numerical analysis will be 



S(T < T rh 
S(T rh ) 



1.01 



(17) 



D. Analytical versus numerical 



In Fig. [2] we show the numerical (solid) and analytical 
(dashed) contour lines of A as function of r and p % x /Tf, 
The analytical lines are produced with g(T) = 106.75 
for simplicity while the numerical curves are calculated 
with the exact g(T) of the SM [TJ. As expected from 
Eq. (fT5")) . once T is high enough to guarantee the initial 
radiation-dominated epoch, A is basically independent 
of it for p x /Tf and <?j fixed. 

We also determine numerically the parameter space 
violating the BBN bound ([TB")) . This region corresponds 
to the filled area of Fig. 2, whose border curve T r h — 
Tbbn reaches the line A = 1.01 asymptotically. This 
constraint on r is by one order of magnitude stronger 
than the one obtained analytically by Eq. (|T4"|). 

The numerical analysis highlights the existence of 
small entropy dilution even in absence of a A-dominated 
epoch. Indeed below the dotted-dotted-dashed line of 
Fig. [2] the X field decays before px reaches pn but pro- 
duces anyway a dilution A < 2.8. Approaching this 
regime corresponds to considering t e ~ r for which the 




6 9 12 15 
LogmlVGeV" 1 ] 



FIG. 2. Contours plot of A = 1.01, 10, 10 2 , 10 5 , 10 8 , 10 10 (la- 
bels on the right) as function of Px/T? and r determined an- 
alytically (dashed line) and numerically (solid line) using SM 
degrees of freedom. Below the dotted-dotted-dashed curve 
the X decay injects entropy even without any X-dominated 
epoch. The filled region is excluded by BBN. 



analytic approximation fails, as we already have stated in 
Section TlIB 1 As a consequence, the constraint A < 1.01 
evaluated analytically allows for values of p x /Tf that 
actually are excluded by two order of magnitude. In- 
stead for A = 10,10 2 ,10 5 ,10 8 ,10 10 the analytical and 
numerical curves are in good agreement and the differ- 
ence appearing at large r is due to the simplification 
g(t) = 106.75 in the analytic estimates. 

At this point one might be interested in knowing by 
what factors one should correct the analytic results in or- 
der to reproduce the numerical outcomes. These factors 
can be extracted from Fig. [3] where the absolute percent- 
age errors 



ry irrinum) /rr ,(ana) , i 

Z l = Mrh / 1 rh ~ L \ > 

Z 2 = \{p x /Tf)^/{p x /Tf)^ m ) 



II 



are plotted as function of A (at r = 10 14 GeV _1 but 
Z\ and Zi are essentially independent of r). It turns 
out that the analytic approach overestimates T r h by a 
factor j3 ~ 2.5, typically. Moreover for a given value 
of p l x /Tf. it underestimates the entropy dilution when 
A < 6 and its error remains below ~ 10% when A > 
3. In particular, the error in estimating the dilution by 
Eq. (fl~3"| or (fT5|) cannot be attributed univocally to a 
mismatch in T r h [which might be due to a wrong choice 
of the prefactor in Eq. (fT2")l ] because in such a case the 
equality Z\ = Z2 should always arise (T e is evaluated 
very precisely in both approaches for A> 1). 
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FIG. 3. The corrective factors Z\ 



jn(num) irp(a,na) 
1 rh 



/T^'-M (solid 



line) andZ 2 = \{f? K l T l) {ana) I {A /Tf) {num) -1\ (dashed line) 
for a fixed t as a function of A. The argument in the modulus 
of Z\ (Z2) is negative for A > 1 (positive for A > 6). 



FIG. 4. Contour plots of the bounds A and A(Tg = 10 5 GeV) 
for Y^/F^ p =10,10 2 ,10 5 ,10 6 are marked by solid and dot- 
ted lines, respectively. The region excluded by each bound 
stands above the corresponding curve. 



III. BARYOGENESIS BOUNDS 

The formation of primordial elements requires Ys(t) 
at BBN times to be compatible with the measures (JTJ 
and ([2"]). In the previous section we used these bounds to 
constrain the history of S(t) by assuming iVs(t) constant 
at T < Tqbn- However, further constraints on the X 
decay can be inferred by knowing the evolution of Ng (t) 
before BBN, as we explain now. 

Let us consider a given baryogenesis mechanism that 
at the temperature Tb generates a B asymmetry Yg ax 
at best. Depending on the relative times at which the 
B asymmetry and entropy dilution are produced, three 
different cases are possible: 

1. the entropy injection occurs exclusively after Tb 
and thus the condition A < Y B nax /Yg Xp is required; 

2. the entropy injection happens exclusively before Tb 
(i.e., T r h > Tb) and then any A is allowed; 

3. the period of entropy injection encloses Tb and 
hence the part of entropy produced after Tb must 
be smaller than Yg MX /Yg Xp ^ 

A(T B ) = S(T Q ^ T ; h) < Y B nax /Yg Xp . (18) 



To avoid confusion in the definitions, we stress that A is the 
dilution due to the whole entropy produced by the X decay, while 
A is only that portion of dilution occurring exclusively after Tg . 



Notice that the parameter space fulfilling Eq. (fT5)) in- 
cludes the region contemplated by the possibilities 1 and 
2. Therefore the parameter region where a given baryo- 
genesis mechanism is compatible with the X decay is 
determined by the bound (|18|). 

As example, in Fig. 4 we present the A constraints 
(dotted lines) for some baryogenesis mechanisms gener- 
ating Y B nax /Yg Xp = 10, 10 2 , 10 5 , 10 6 at T B = 10 5 (the fig- 
ure also reports the numerical lines A = 10, 10 2 , 10 5 , 10 6 
in solid) . The parameter space excluded by each A bound 
stands on the right of the corresponding line 0. 

The curves A and A overlap asymptotically because at 
large r the entropy is injected mostly after Tb and hence 
we have S(T B ) ~ S(Ti) in Eqs. © and (IS]). Thus, in this 
regime the parameter space excluded by Eq. (|18|) can be 
derived analytically by the condition A < Y B nax /Yg Xp . 

Instead, in the region where the curves A are vertical 
there is no bound on A but only on r. In such a region the 
A condition is much less stringent than naively requiring 
T r h > Tb (which actually corresponds to A ~ 1). The 
reason is that big values of Y B nax /Yg Xp allow for large 
dilutions of the produced B asymmetry so that a region 
with T r h < Tb is permitted. This is quantified in Table 
I where for various values of Y™ ax jY^ p we report the 



7 In realistic cases Tg stands for the temperature interval during 
which the B asymmetry is created. Here we work in the limit 
that the entropy produced during this interval is negligible. On 
the contrary, the vertical part of the A curves is stumped while 
its oblique section holds sharp. 



(i 



A < 


10 


10 2 


10 3 


10 4 


10 5 


10 6 


tB/T > 


2X10" 1 


3xl0~ 2 


5xl0~ :i 


7xl0~ 4 


lxl0~ 4 


2xl0~ 5 


T B /T rh < 


3.4 


5.5 


8.8 


14 


21 


33 



TABLE I. Numerical evolution of the bound on r (T r h) as 
function of A and ts (?s) in the regime of very large A. 



bounds on ts /t and Tb /T r h valid in the regime of very 
large A, that is, when the curve A is vertical. 

The interesting issue is that the values in Table I are 
basically independent of variation on tg and Tg. More- 
over the shape of the curves of the bounds A is univer- 
sal, as we see in Fig. 4 where every dotted line could be 
obtained by the translation of one of the others. This 
implies that by our results anyone can take some long- 
lived field and baryogenesis mechanism and determine 
their compatibility without implementing further numer- 
ical analyses. To do it, one needs only to locate the 
curve ([15]) in the plan T-p l x /Tf for the desired frame- 
work. This can be achieved as follows. 

First of all, one has to obtain the horizontal position 
of the A curve. This is provided by the constraint on r 
in presence of very large A (namely, the position of the 
vertical part of the A curve). It can be found by Eq. (fT4"]) 
relating r to the analytic estimate of T r h which in turn 
is constrained by 



> -T B (for large A) 



(19) 



In this expression 7 is the opportune value in the second 
row of Table I connecting Tg to the numerical estimate 
of T r h, while /3 ~ 2.5 is the corrective factor at large A 
determined in Section Hi Dl Subsequently, the horizontal 
position of the A curve is given by 



r < 7.7 x 10 17 GeV- 



(for large A) . (20) 



Afterwards, one needs to fix the vertical position of the 
A curve. This is furnished by Fig. 2 because for r — > 00 
we have 



A 



A < Y^ Qa 



lyexp 



(for large r) 



(21) 



In conclusion, the compatibility of late-time entropy 
injection with successful baryogenesis can be calculated 
by using simple arithmetic. Some concrete examples will 
clarify how this procedure can be applied in realistic sce- 
narios. 



IV. EXPLICIT APPLICATIONS 

The procedure we have just presented is model inde- 
pendent except but the assumptions: 

i) the X decay follows Eq. ([3]), does not induce a B 
asymmetry and its products thermalize fast. 



Moreover it is opportune to have: 

ii) the BAU production is much faster than the en- 
tropy injection. 

Therefore, the described procedure allows for investi- 
gating the compatibility between any specific baryogen- 
esis mechanism and any particular X field satisfying the 
conditions i and (possibly) ii. As explicit applications, 
we will analyze some concrete mechanisms embedded in 
the MSSM whose g(T) is calculated taking the scalar 
[fermion] content at 0(1) [0(0.1)] TeV for definiteness. 



A. Some baryogenesis mechanisms 

Here we apply our results to electroweak baryogene- 
sis and thermal (resonant and no 11 resonant) leptogenesis 
in presence of a generic X field. For each baryogene- 
sis framework we first review the estimate of Y B nax 



land 

then we calculate the parameter region where the X field 
does not destroy the BAU. The results will be summa- 
rized in Fig. 5, in which also the BBN bound and some 
A contour (dotted dashed) curves for the MSSM are re- 
ported. 



1. Electroweak Baryogenesis 

In electroweak baryogenesis the observed BAU is pro- 
duced during the electroweak symmetry breaking [18|. 
In this scenario the departure from thermal equilibrium 
is achieved by a first order electroweak phase transition 
(EWPT) . In such a case the transition proceeds via nucle- 
ation of bubbles containing the electroweak broken phase 
and the movement of the bubble breaks locally the ther- 
mal equilibrium conditions. Then, just in front of the 
expanding bubbles, C and CP violating interactions gen- 
erate a left-handed asymmetry that SU (2) j, sphalerons 
transform to a B asymmetry entering the bubble. How- 
ever, the formation of the B asymmetry is not enough: 
it is also necessary to preserve it till today. This occurs 
if the first order EWPT is strong @. 

In the SM the EWPT is not strong and the CP violat- 
ing sources are too small to form enough B asymmetry 
[18| . These problems are overcome in some extensions of 



8 We take the estimates of the literature assuming standard cos- 
mology but in principle the evaluation of Yg ax could change 
when the BAU is produced in the X-dominated era. However, 
for our concrete baryogenesis examples, Refs. |l6t \vf\ find small 
modifications of YJ^ ax and Tb so that we can use the standard 
results. 

9 The EWPT is strong when SU(2)l sphalerons are out-of- 
equilibrium inside the bubbles. Unless of subtle circumstances 
Il7l. Il9t not considered here, this happens for v(Tg)/Tg > 0.7 
where v = v(T = 0) = 174 GeV and v(Tb) is the vacuum expec- 
tation value of the (SM-like) Higgs at the temperature Tg when 
the phase transition begins Il8ll - 
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FIG. 5. Numerical contours (dashed dotted lines) of A = 
1.01,10, 10 s , 10 s (labels on the right) as function of p l x /T? 
and t for MSSM degrees of freedom. Regions excluded by 
BBN, electroweak baryogenesis (EWBG), thermal leptogene- 
sis (TLP) and resonant thermal leptogenesis (RTLP) are filled 
as reported in the legend. 



the SM. Some non-supersymmetric scenarios have been 
considered but at present the analyses of them are 
not so developed to provide precise estimates of the maxi- 
mal B asymmetry these models can produce. Instead ac- 
curate predictions exist for supersymmetric extensions. 

The MSSM can reproduce the observed BAU if the 
gaugino-Higgsino sector is at the electroweak scale 
|23 |. the SM-like Higgs and right-handed stop are light 
(m h < 127 GeV, m iR < 120 GeV) and the left-handed 
stop is heavy (mj > 6.5 TeV) [26j . Moreover, the B 
asymmetry formation typically starts at Tg ~ 125 GeV 
and ends after a few GeV [26|. Finally, concerning Yg iax , 
different treatments of CP- violating sources and flavor 
effects exist (23l - |25l | and they lead to B asymmetries that 
may differ by almost one order of magnitude. Choosing 
the intermediate result we can consider Y^ ax ~ 10 Yg Xp 
[2IEH]. 

In conclusion, the parameter space where electroweak 
baryogenesis in the MSSM is compatible with the X de- 
cay is given by the constraint 

A(T B ~ 125 GeV) < 10 , (22) 

in which the condition ii is fulfilled because the B asym- 
metry is produced in a narrow temperature interval 
0(1 GeV). 

The numerical solution of Eq. ([22]) gives the area la- 
belled EWBG in Fig. 5. Its border curve could be cal- 
culated by the analytic procedure expressed in Section 
IIIII In fact, the position of the vertical part of the curve 



is given by Eq. flU) yielding r < 5 x 10 13 GeV -1 for 
9rh = g(T r h = 92 GeV) ~ 100 and 7 = 3.4, while at 
much larger r the line must coincide with the bound 
A < 10 because of Eq. (f2~T|) . Instead, in the interme- 
diate regime of r the curve is deduced by the universal 
shape of the A lines of Fig. 4. 

Before concluding, some words about electroweak 
baryogenesis in non-minimal supersymmetric extensions 
are in order. The estimate Yjp ax Pzil0Yg Xp holds roughly 
correct also in extensions of the MSSM where the extra 
content couples weakly to the Higgs sector and no new 
large sources of CP violation are introduced. This is not 
the case if an extra singlet is added to the MSSM as it is 
performed in Ref . [27j ■ In such a extension Y^ ax can be 
enhanced approximatively by a factor 5. Moreover, other 
modifications of the minimal setup can further enhance 
Y]^ ax , as for instance in the extension called Beyond- 
the-MSSM where Y^ ax may increase by a further or- 
der of magnitude [28j. Here we do not explicitly repeat 
the analysis for all these variations but their compati- 
bility with a late-time entropy injection could be easily 
figured out algebraically as we have just performed for 
the MSSM. 



2. Thermal Leptogenesis 

Leptogenesis is one of the attractive scenarios to ex- 
plain the BAU [29|, |3(j. Its central idea is that in 
a S-symmetric but Lepton(L)-asymmetric Universe it 
is possible to produce the BAU thanks to electroweak 
sphalerons that equilibrate the B+L asymmetry with- 
out changing B—L. In thermal leptogenesis the required 
initial L asymmetry can be obtained for instance in the 
Type-I seesaw framework (3l| where one singlet right- 
handed neutrino TV per family is added to the SM 0- 
These fields have the interactions 

C D ^(M N ) u N t N t + yijNiljiTiH* + h.c. , (23) 

in which I (if) is the SM lepton (Higgs) and indeces run 
over families. They not only give rise to a net L asymme- 
try in the early Universe by out-of-equilibrium decays but 
also generate sub-eV neutrino masses via the canonical 
seesaw mechanism as required by the neutrino oscillation 
data [Hj]. This mass turns out to be m v w \y\ 2 v 2 /Mjy 
yielding M N ~ 10 14 GeV for \y\ ~ 1 and m v ~ 0.1 cV. 

Assuming a normal hierarchy in the heavy neutrino 
sector, the CP asymmetry of the decay of the lightest 
right-handed neutrino N\ is given by [34| 

|£/| = l^v / ^^ sin ' 5 ' ( 24 ) 

where Am^ tm is the atmospheric mass scale of light neu- 
trinos [13 and 8 is the effective CP violating phase. 



Our conclusions do not change for Type-II and Type-Ill V't'M . 



s 



The L asymmetry, which is mostly produced by the 
decay of N\, is given by = e I Y^ 1 W where the ther- 
mal wash out W takes into account the effect of inter- 
actions reducing the created L asymmetry, as for in- 
stance IH -H* IH* . Afterwards, electroweak sphalerons 
tend to equilibrate the B+L asymmetry and so they con- 
vert around half part of Yl into Yg. Therefore, in order 
to determine the maximal B asymmetry YJ^ ax , one sets 
sin 8 = 1 and obtains 1 341 



1 B 

v exp 
1 B 



4 x 10 4 



/ m Ni > 


I V /Am atm 


( Y Nl W \ 


I 10 14 GcV J 


0.05eV 


{7x 10- 4 J 



(25) 

Since commonly the framework is supposed to be embed- 
ded in a grand unified theory and moreover the Universe 
typically never reached temperatures high enough to 
thermally generate neutrinos beyond the GUT scale, we 
can consider 10 16 > M^GeV -1 > fO 9 , where the lower 
bound comes from requiring enough BAU in Eq. (|2"5"1) . 

Now we can analyze the B-asymmetry dilution due to 
X decay EE Taking T B « M Nl , the case M Nl ~ 10 9 GeV 
gives Yg lax /Yz p = 0(1) and yields then 



A(T S ~ 10 9 GeV) < 0(1) 



(26) 



Eq. ((261) permits any A for r < 10 -1 GeV -1 but does 
not provide a sharp bound on r because the condition ii 
is not guaranteed. Similarly for Mn 1 ~ 10 16 GeV, which 
yields Yg 1 ^ /Yg Xp = O(10 6 ), we have the constraint 



A(Tb ~ 10 ib GeV) < C(10 b ) 



(27) 



that allows any dilution A < 10 6 independently of r. 

The union of the regions fulfilling (|26|) and (|27|) gives 
the parameter space where thermal leptogenesis is com- 
patible with the X decay. The numerical result is shown 
in Fig. 5 (green area labelled TLP) and could be repro- 
duced algebraically by the procedure explained in Section 
IIIII and already applied for electroweak baryogenesis. 



3. Resonant thermal leptogenesis 



The CP asymmetry in Eq. (|24[) requires very heavy 
additional neutrinos to produce the BAU. For this rea- 
son thermal leptogenesis with normal neutrino hierarchy 
provides no experimental evidence at current achievable 
energies. This unappealing feature is avoided in resonant 
thermal leptogenesis where extra neutrinos at the elec- 
troweak scale may yield BAU and detectable signatures 
at the same breath (35j. 

The key idea of resonant leptogenesis (36l . [37j is that 
the self-energy effects dominate the leptonic asymmetries 



11 An analysis focused on solving the gravitino problem in a par- 
ticular framework of thermal leptogenesis and late-time decay is 



presented in Ref. 



when the mass splitting between the right-handed neutri- 
nos N is much less than their masses. If the mass splitting 
between these fields is comparable to their decay widths, 
the CP asymmetry gets enhanced resonantly. 

Focusing on type-I leptogenesis (our conclusions would 
not change for type II and III) , let us consider two singlet 
Majorana fields N\ and iV 2 of masses Mi and M2. Their 
Yukawa interaction hijNiljH allows for decays to SM 
lepton I and Higgs doublet H whose CP asymmetry is 
given by [H] 



ei = 



Im(/it/i)2 2 



{Ml 



M%)MiM 2 



8 7 r(/it/ l ) 11 {Ml - MlY + (M 2 r 2 - Miro 



(28) 



Hence, assuming Mi ~ M 2 and Mi — M 2 ~ T\ — T 2 we 
can achieve e\ ~ 0{1) so that a large L asymmetry Yl 
can be produced even if the initial N\ and N 2 abundance 
Yv is small. Finally, SU{2)l sphalerons can convert Y^ 
to Yb if the lepton asymmetry is produced enough before 
T ~ 130 GeV when sphalerons decouple (39j. Typically, 
this happens if Mi > 250 GeV (37j . In such a case it 
turns out Yb — 0.5 x eiYjyW, where it is reasonable to 
consider Y N W = O(10" 4 ) (H| yielding Yg lax = O(I0~ 4 ) 
for ei - O(l). 

In conclusion, in resonant leptogenesis it turns out to 
be Y B nax w O{10 6 )Yg Xp . Moreover, in order to favor the 
compatibility with the late-time decay, Tb has to be as 
low as possible. Since it is Tg w Mi, the most favorable 
choice is Mi ~ 250 GeV. Subsequently, resonant leptoge- 
nesis is compatible with the X decay when the constraint 



A(T B ~ 250 GeV) < O(I0 b ) 



(29) 



is fulfilled. Observe that this constraint can be considered 
as a sharp bound since the condition ii is satisfied due 
to the short temperature interval during which the B 
asymmetry is generated (250^T/GeV^ 130). 

By solving Eq. (|2"9")l numerically, one finds that the ex- 
cluded region is the area labelled RLPT shown in Fig. 5. 
As already checked for the previous baryogenesis mech- 
anisms, the excluded region could be easily determined 
algebraically by the procedure of Section IIIII Indeed its 
vertical border is due to the bound t^Ix 10 15 GeV -1 
coming from Eq. fl20|) with g rh = g(T rh = 18 GeV) ~ 90 
and 7 = 33, its lower part corresponds to the constraint 
A ~ A < 10 6 shown in Fig. 2, and the curved part is 
deduced by copying the shape of the A lines of Fig. 4. 



B. 



X field as a modulus 



In the examples we have just considered the field pro- 
ducing the late-time entropy is not specified. Now we 
analyze the particular case of the X field being a modu- 
lus. Let us assume it to be gravitationally coupled to the 
SM thermal bath via a dimension-five operator so that 
its lifetime is given by 



2tt 



(30) 
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FIG. 6. Implications of the analysis summarized in Fig. 5 for 
decay of muduli with coupling ax = 1. Numerical contours 
(dotted dashed curves) of A = 1.01, 10, 10 5 , 10 8 (labels on the 
right and top) as function of Xi and Mx- Regions excluded 
by BBN, electroweak baryogenesis (EWBG), thermal lepto- 
genesis (TLP) and resonant thermal leptogenesis (RTLP) are 
filled as reported in the legend. 



where ax spans from 0(l)toO(lCT 2 ) depending on the 
non-renormalizable coupling (40| . 

In an expanding Universe, when the Hubble scale is 
roughly equal to the modulus mass, H ss Mx, the mod- 
ulus enters into the oscillating regime. By taking the 
initial amplitude of oscillations of X to be Xi , the energy 
density p x can be expressed as 



fix 



M 2 x Xf 



(31) 



Moreover, when the field starts oscillating, the Universe 
is radiation dominated and its temperature T, is 



T 



90 m 2 P M% 
k 2 9i 



(32) 



Subsequently, we can easily convert the model indepen- 
dent bounds of the previous sections to the case of X as 
a modulus. For instance, the general bounds on r and 
p x /Tf presented in Fig. 5 can be re-expressed as con- 
straints on Xi and Mx, as shown in Fig. 6. 

Using Eqs. (|30 |) -([32 | . the sudden decay approximations 
©, (JT3J) and HU) become 



ji4 



10 XfM x 
9ir 2 gi m & p 

10 a 2 x M x 
7r 4 3r/ l mp 



(33) 
(34) 



y 9 gt s/aHmpMx 

and the BBN bound T rh > 4 McV is converted to 

M x > 7 x 10 4 (l/a x ) 1/3 GcV (analytic) , (36) 

which is slightly weaker than the numerical one. How- 
ever, the BBN bound is not the strongest constraint in 
the regime of large oscillations. In such a regime, to avoid 
the BAU wash out, we need M x ^/aT > O(10 13 ) GeV for 
thermal leptogenesis, Mx-&ctx > O(10 8 ) GeV for elec- 
troweak baryogenesis and Mx^Jctx > O(10 7 ) GeV for 
resonant leptogenesis. 

V. CONCLUSIONS 

There exist many theories beyond the standard model 
of particle physics that predict cosmologically long-lived 
fields. Even though the initial energy density of these 
fields may be many orders of magnitude less than the 
one of radiation, at later times they can dominate the 
expansion of the Universe and dump sizeable amount of 
entropy during their decay. Subsequently these fields di- 
lute the pre-processed B asymmetry and may alter the 
primordial element abundances. 

In order to constrain the above theories, in this paper 
we analyzed the entropy and the subsequent re-reheating 
that a generic long-lived field X produces. We solved the 
decay equations numerically and we checked their con- 
sistency with the analytic sudden-decay approximation. 
The result can be parameterized by lifetime and initial 
energy density of the X field. We found that for small 
entropy dilution A the analytic approximation badly fails 
because in this regime there is no clear A-dominated 
epoch. Instead for A > 3 the analytic and numerical 
calculations disagree by less than 10% in the entropy es- 
timates but by a factor ^2.5 in the evaluation of T r h, 
the temperature at the end of the dilution. However, 
the errors of the analytic approach can be overcome by 
applying some corrective factors which, among other ef- 
fects, strength slightly the BBN upper bound on the X 
lifetime. 

The revised control on entropy dilution and re- 
reheating temperature allowed us for investigating the 
effect of the X decay on the B asymmetry. The analysis 
was carried out by assuming the mechanism producing 
the entropy injection to be independent of the one gen- 
erating the BAU. The outcome is a set of conditions on 
the lifetime and initial energy density of X as function of 
when and in what amount the BAU is produced. These 
conditions are quite generic so that they can be easily ap- 
plied to determine the compatibility of baryogenesis and 
late-time decay described by a wide class of scenarios. 

As illustrative examples, we applied our results to 
some concrete models. It turned out that success- 
ful MSSM electroweak baryogenesis, thermal leptoge- 
nesis and resonant leptogenesis put strong bounds on 
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the initial abundances of X fields with lifetime r > 
5 x 10 13 ,10" 2 ,10 15 GeV -1 , respectively. Instead, for 
smaller values of r the abundance is not constrained (see 
Fig. 5) . Furthermore, if we take X being modulus as ex- 
plicit example of decaying field, the above constraints be- 
come bounds on the initial oscillation amplitude X$ , mass 
Mx and gravitationally-mediated coupling ax to the vis- 
ible sector. In particular, very large Xi implies excessive 
dilution of the B asymmetry produced by MSSM elec- 
troweak baryogenesis, thermal leptogenesis and thermal 
resonant leptogenesis for M x ^fax > 10 8 , 10 13 , 10 7 GeV, 
respectively (see Fig. 6). 
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APPENDIX: RELAXING THE BBN BOUND 

Apart from some special exceptions 0, the observed 
primordial element abundances cannot be explained 
when sizeable entropy injections occur during the BBN 
epoch. On the contrary, the success of the standard BBN 
model is not spoiled for Q 



where the A-decay hadronic channels are supposed sup- 
pressed. 

In our numerical analysis we implement the reheating 
temperature as 



S{T < T rh ) 
S(T rh ) 



(38) 



T rh > T BBN = 4 McV 



(37) 



with £ = 0.01. Subsequently T rh is the temperature of 
the thermal bath when the 99% of the entropy due to the 
X decay has been injected. 

Applying the definition (|3"B"|) in the constraint (|3"T|) is 
consistent with BBN. Indeed, even in the pathological 
combination of baryogenesis and entropy mechanisms 
leading to T rh = T BBN and Y B (T BBN ) ~ 8.4 x lO" 11 
[the minimal value allowed by Eqs. |T]) and ([2])], the re- 
maining 1% of entropy injection after Tbbn dilutes Y B 
by an amount that is practically negligible for the exper- 
imental constraints. 

However, the convention £ = 0.01 is often too conser- 
vative. For instance, in the case with T r h = Tbbn and 
Yb(T BB n) — 9.2 x 10~ n , the maximal allowed entropy 
dilution after T BB n is the ratio between the upper and 
lower bounds of Y^ xp . In such a case, by taking the al- 
lowed extrema of ([1]) and @, the BBN bound can be 
implemented by (f3"T)) where T r h is redefined through (|3"5)) 
with £ = 0.10. On the other hand, if possible system- 
atic errors due to priors were taken into account in the 
WMAP analysis, the exp erimental bound © would be 
comparable with (JXJ) [12| and the BBN bound could be 
further relaxed by redefining T r h by (|3"8|) with £ = 0.27. 

As a consequence, these different choices of £ slightly 
relax the BBN bound shown in the figures of the paper. 
For instance, the convention £ = 0.27 weakens the BBN 
constraint on p x /Tf (on r) by a factor ~10 (~3). In any 
case for the possible values 0.01 < £ < 0.27 the numerical 
analysis still provides a BBN bound that is stronger than 
the one obtained analytically by Eq. (|14|). 
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